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Abstract 

A model of sound generated in a high subsonic 
(Mach 0.9) circular jet is solved numerically in 
cylindrical coordinates for nonaxisymmetric distur- 
bances. The jet is excited by transient mass injec- 
tion by a finite duration pulse via a modulated ring 
source. The nonaxisymmetric solution is computed 
for long times after the initial disturbance has exited 
the computational domain. 

The long time behavior of the jet is dominated 
by vorticity and pressure disturbances generated at 
the nozzle lip and growing as they convect down- 
stream in the jet. These disturbances generate sound 
as they propagate. The primary nonaxisymmetric 
effect that we simulate is that of a flapping mode 
where regions of high and low f pressure alternate on 
opposite sides of the jet. The predominant feature of 
this mode is the appearance of relatively large devi- 
ations of the pressure from the ambient pressure on 
opposite sides of the jet and the convection of these 
regions downstream. 

We illustrate flow field, near field and fax field 
data. Important nonaxisymmetric characteristics of 
the near and flow field disturbances include roughly 
periodic pressure elevations and depressions at op- 
posite values of the azimuthal angle t/u These cor- 
respond to pressure disturbances propagating in the 
axial direction. The azimuthal velocity exhibits a 
sinusoidal dependence on with similar roughly pe- 
riodic disturbances. For every azimuthal angle ‘ip 
the jet radiation peaks about 30° from the jet axis, 
however there is now a pronounced dependence of 
the far field radiation pattern on rj>. 
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1 Introduci ion 

In this paper, the generati m and propagatio i 
of sound in a high subsonic jet excited by nonax- 
isymmetric disturbances is simulated. We consider 
a circular jet exiting from a straight nozzle extend- 
ing to infinity in the upstream cirection. The jet is 
excited by a nonaxisymmetric transient source and 
the long time behavior of the pressure is considered. 
We consider the full nonlinear Euler equations in 
conservation form modified so that a spreading je 
profile, obtained from experimental measurements, 
is a solution in the absence of any disturbances. 

The primary objectives are (i) o characterize non 
axisymmetric aspects of the flovv field and near field 
and (ii) to characterize nonaxisymmetric effects in 
the far field. The jet is initially excited by a spatially 
and temporally localized source of transient mass in 
jection. In contrast to previous work we consider 
modulated ring source for the excitation, thus gen 
erating a nonaxisymmetric disturbance field. Thi 
leads to an initial acoustic distu bance which prop 
agates through the jet. As a result of the excitation 
instability waves are generated in the jet. Thes* 
waves grow and then decay as t hey convect down 
stream, generating sound in the process. This phe 
nomenon occurs over time scale" much longer thai 
that of the excitation pulse. T\ e emphasis in thu 
paper is on the long time response of the jet, after 
the wave due to the excitation pulse has exited tht 
domain of interest. Thus, instabi ity wave generated 
sound in the jet is simulated. T he simulation doer 
not directly account for sound generated by smali 
scale turbulent sources in the jet 

In previous work we have co isidered jet acous 
tics for circular jets subject to axisymmetric dis 
turbances. In particular, both liuear and nonlineai 
responses have been considered L * 2> 3 In these ref 
erences the jet is excited by a ;ource of transient 
mass injection and the development of instabilities 
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together wit.li associated sound generation is com- 
puted. 1 ho results in these references demonstrate 
the properties of instability wave generated sound. 
A succession of disturbances is generated at the noz- 
zle lip. These disturbances initially grow as they 
convcct in the jet. As they grow they generate pres- 
sure disturbances which propagate into the far field 
as sound. The far field radiation peaks at about 30° 
from the jet axis in agreement with measurements 
and observations. 

We have also considered two dimensional Carte- 
sian jets. 4, 5t 6 * 7 In two dimensions we have com- 
puted both the sound generated from instability 
waves and the interaction of this sound with a nearby 
flexible structure . These computa- 

tions spanned a range of jet velocities up to and 
including supersonic jets. The two dimensional re- 
sults for high subsonic jets (the parameter range 
considered here) indicated that the long time be- 
havior of the excited jet is dominated by a nearly 
periodic acoustic response at a particular fundamen- 
tal frequency together with harmonics. 6 This radi- 
ation is beamed at approximately 30° from the jet 
axis. The unsteady flow field is dominated by distur- 
bances generated at the nozzle lip and propagating 
along the jet. In contrast to the nearly periodic far 
field the flow field exhibited an aperiodic, apparently 
chaotic behavior. 

The exact sources of jet noise have been 
identified from the basic equations of fluid 
dynamics. 8, 9 > 10, 11 * 12 Computation of the exact 
sources would require solution of the Navier-Stokes 
equations in the jet accompanied by an appropri- 
ate turbulence model. Such computations are diffi- 
cult to conduct, particularly in 3 dimensions. Both 
small scale and large scale disturbances in the jet 
can act as sources of sound. In this paper we con- 
sider only sound generated from large scale insta- 
bility waves. These instability waves are generated 
within the transiently excited jet by the mechanism 
described above. 

Instability waves or large scale structures which 
can act as sources of sound 
in a jet have been observed in experiments 1, 13 ’ 14 
and studied by analytical 15, 16, 17, 18, 19, 20 - 21 and 
numerical 1, l ’ 3 * 22 methods. While there have been 
a limited number of 3 dimensional computations 
23, 24 > the cost of the 3D computations has to 
date precluded a detailed study of the role of non- 
axisymmetric large scale structures as sources of jet 
noise and has generally forced restrictive assump- 
tions such as linearity. On the other hand, many 
of the large scale structures observed in experi- 
ments are intrinsically nonaxisyinmetric, e.g., flap- 


ping modes where regions of large pressure alternate 
on different sides of the jet, and thus a full under- 
standing of jet acoustics requires an understanding 
of nonaxisymmctric modes in the jet and the result- 
ing acoustic radiation. 

In our model we solve a modified version of the full 
nonlinear Euler equations in cylindrical geometry. 
The equations are modified to account for a spread- 
ing jet as described below. As a result, the invis- 
cid, instability wave sources of sound are computed 
directly together with the resulting sound genera- 
tion. We specifically include a straight nozzle from 
which the jet exits. The flow field and near field dis- 
turbances are computed together with the far field 
sound radiation. 

2 Model and Numerical Method 

The computational domain is shown in Figure 
1. We consider a jet exiting from a straight, circu- 
lar pipe of diameter D , where D = 4 inches. All 
boundaries are artificial boundaries at which radia- 
tion boundary conditions are employed. 1, 27 The jet 
is excited by a transient source of mass injection lo- 
cated approximately ID downstream from the noz- 
zle exit. This source is both spatially and temporally 
localized, leading to a broad band excitation. Un- 
steady pressure, density and velocity are computed 
both interior and exterior to the pipe. This leads to 
the generation of a train of instability waves which 
propagate along the jet, decaying beyond the poten- 
tial core of the jet and generating sound. 

The Euler computations employ a (2-4) version of 
the MacCormack scheme. 25 Further details on the 
numerical scheme are given in the references. 4, 5> 6 
The spatial fourth order accuracy is necessary in or- 
der to follow the acoustic waves over large distances 
without undue numerical dispersion. This is partic- 
ularly important given the coarse grids necessitated 
by nonaxisymmetric computations. 

The Euler equations are solved in conservation 
form in cylindrical coordinates z,r, where z de- 
notes the axial distance, r the radial distance and 
rp is the azimuthal angle. The dependent variable is 
the vector 

q = (p,pu,pv,pw,E) T , 

where p is the density, u, v and w are the axial, radial 
and azimuthal velocity components respectively, and 
E is the total energy per unit volume. The pressure, 
p, is obtained from the equation of state. 

The Euler equations are modified in the jet do- 
main to account for the jet flow. The jet exits from 
a nozzle of diameter D located at ; = 0 and the 
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solution is computed both within and exterior to 
the nozzle. The Euler equations are modified to 
account for two different non- homogeneous forcing 
terms. 6 One term serves as an excitation pulse to 
excite the jet. It corresponds to a localized source 
of mass injection. In axisymmetric computations we 
take a point source centered on the axis at z s . In the 
nonaxisymmetric computations we consider a ring 
source of radius r s , where r 9 ~ . 3D . For each ip 
the source is a peaked Gaussian centered at r Sl z 3 . 
In order to obtain a nonaxisymmetric excitation the 
strength of the source is modulated by the function 
1 + e cos ip, where 6 = 0.2. The temporal dependence 
of the source is described in the references, e.g., 4 . 
The source is localized in time leading to broad band 
excitation. 

The second forcing term is designed so that in the 
absence of the excitation pulse the solution to the 
Euler equations would be a stationary profile cor- 
responding to a spreading jet. In the present com- 
putation, mean density and pressure are assumed 
constant, as are the mean radial and azimuthal ve- 
locities. The mean axial velocity U is a taken from 
experimental measurements and models an isother- 
mal jet with exit Mach number 0.9. 26 

The inclusion of this term separates the compu- 
tation of the disturbance, in particular the result- 
ing instability waves, from the computation of the 
mean flow (i.e. the spreading jet) . Thus, the result- 
ing system of equations allows for the simulation of 
instability waves and the resulting sound generation, 
together with the propagation of acoustic waves in 
the jet flow field, without requiring the computation 
of the spreading jet itself. Although this is a simpli- 
fied model, it captures many of the observed features 
of instability wave generated jet sound and permits 
high resolution computation of the coupling of jet 
noise with, the flexible panels and the resulting radi- 
ation from the panels. In particular, the model al- 
lows for computation of many of the natural sources 
of jet noise (the instability waves) together with the 
sound radiated by these sources. 

3 Numerical Results 

The primary objective of the computations 
is to examine flapping modes, in which regions of 
high pressure in the jet alternate at opposite az- 
imuthal angles. We consider the unsteady pressure 
p — p{t , r, z, ip) - ^ where p Q 0 is the ambient pres- 
sure. In response to the excitation source an ini- 
tial acoustic wave emanates from the source loca- 
tion through the jet. In addition, a large slowly 
moving structure, an instability wave, propagates 


downstream along the jet axis We consider the e 
transient effects. We concentrate on the long tine 
behavior of the jet which is dom inated by by a shed- 
ding of disturbances from the j.ozzle lip, the prop- 
agation of these disturbances ( ownstream and the 
generation of sound from these disturbances as they 
grow and decay along the jet. We note that theie 
are also disturbances that propagate upstream insice 
the pipe. 

In Figure 2 we consider cok r contours for p i n 
the r — z plane sin(^) = 0, i.e , the pressure field 
for the angles ip = 0 and ip = 7r. The coloring s 
according to the spectrum with red indicating large 
positive values of p and violet in dicating laxge nega- 
tive values of p. Note that color are assigned based 
on a logarithmic scale. These contours are shown in 
Figure 2 at 6 different times, t - , . . £ 6 , where the 
nondimensional time t is define d by £ = iU e /D , £ 
is dimensional time and U e is i le jet exit velocity. 
The times illustrated in Figure 2 are t- L = 28. 0£, 
t 2 = 29.19, t 3 = 30.83, t A = 31.93, t 5 = 33.02, 
£ 6 = 33.56. 

Note that these figures show e lly a portion of th 3 
larger computational domain, /hich extends 30 D 
both upstream and downstream of the nozzle exic 
and radially. In interpreting the ;e figures, we recall 
that the nozzle exit is at z = 0 a id the initial sourc 3 
was a ring source centered appn dmately ID down 
stream of the nozzle exit and of pproximately 0.51 f 
in radius. 

For each time there is am array of acousti : 
waves propagating primarily up scream and emanat 
ing from disturbances near the nozzle exit. In addi 
tion, large radiation, predominantly yellow in color 
can be seen emanating at mid angles from the jet 
Note that the mid angle radiatic n is of greater spa 
tial extent than the upstream radiation, indicating 
a low frequency, relatively broac band spectrum ai 
mid angles and a higher frequency, sharper spectrum 
radiation field at upstream angles. In addition, th< 
propagation of disturbances up the nozzle can b( 
clearly seen. The flapping mod^ is best visualized 
by the alternating regions of larg radiation (yellow 
emanating from disturbances in the jet. At time: 
t\ and £2 this radiation is most pronounced in the 
lower half of the figure (ip — 1 r). whereas for timee 
£3 — £6 the predominant near fir id pressure distur 
bance is in the upper portion of the figure (ip — 0) 
As time evolves these pressure disturbances propa 
gate into the far field and weak n due to spherica 
decay. Thus, for example, the p = 0 yellow dis- 
turbance seen at approximately ID downstream at 
£1 becomes the yellow-green dist ir bance at approx 
imately 5 D downstream at t 2 wl ile a yellow distur- 
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bance has been generated in the lower half of the 
figure {xp = 7r ) figure at approximately 2D at time 
£3. A similar effect is seen at other angles, although 
the greatest contrast is for the xp = 0/xp = n plane, 
which is also the plane of the largest asymmetry of 
the excitation source. Note that the figure also in- 
dicates that pressure disturbances emanating from 
the shear layer become coupled with the rigid pipe 
as they propagate outward. The coupling switches 
from top to bottom as the disturbances in the shear 
layer switch. 

Another perspective of this flapping is shown in 
Figure 3 where pressure in a cross section of con- 
stant z/D is shown at 6 different times, t a , . . . , tf. 
Here, t a = 22.61, t b = 24.26, t c = 24.80, t d = 27.00, 
t € = 28.09 and tf = 33.02. The cross section is 
taken at the axial location z/D — 4.89 just before 
the end of the potential core of the jet. The diameter 
of each cross section is approximately 2D, the col- 
oring is again according to the spectrum, and, due 
to the wide range of values involved, for each time 
the data is internally scaled. For each of cross sec- 
tion xp = 0 corresponds to the right horizontal axis 
while xp — 7r corresponds to the left horizontal axis. 
Thus, at times t a and t b there is high pressure near 
%p = 7r for this particular axial location, while at 
times t e and tf the high pressure occurs near xp = 0 . 
From Figure 2, it can be inferred that there are dis- 
turbances with high pressures alternating between 
ip = 0 and xp = 7r shed from the nozzle lip and prop- 
agating downstream along the jet. It can also be 
seen from Figure 3 that the solution is symmetric 
with respect to reflections across the horizontal axis 
in the cross section, i.e., under the transformation 
rp —xp. The nonaxisymmetric source distribu- 
tion exhibits this symmetry and it appears to be 
maintained by the solution. 

We next examine the flow field azimuthal velocity 
w. We consider w as a function of xp at a particu- 
lar flow field location near the edge of the jet shear 
layer. In Figure 4, w is plotted as a function of ip for 
a range of values of t (after the large instability wave 
has passed this point). The data is taken at approx- 
imately 3,5 diameters downstream and inside the jet 
flow (r/D = 0.7). It can be seen from the figure that 
as a function of t p, xu exhibits a cosine dependence 
(inherited from the initial nonaxisymmetric distur- 
bance) together with ripples that appear at roughly 
equally spaced intervals in time. These ripples rep- 
resent local maxima and minima of w in time and 
are associated with a succession of flow field distur- 
bances generated at the nozzle lip and convecting 
downstream in the jet (these disturbances can also 
be seen in Figure 2). It appears that the angular 


locations of these ripples (i.e., of the maxima and 
minima of disturbances passing a fixed point) un- 
dergoes a slow rotation in ip as can be seen from the 
lines drawn on the figure. 

In Figures 5a and 5b we examine p outside of the 
flow field at a distance 6D from the source center. 
Examination of points at greater distances from the 
source center indicates that this figure qualitatively 
represents far field acoustic pressure. In Figure 5a, 
p is plotted against t for points at 30° from the jet 
axis for xp = 0 and xp = n. Figure 5b is similar 
except that it is for points at 90° from the jet axis. 
Both figures use data only for times after the leading 
acoustic pulse has passed the given points. 

The figures at both points illustrate the effect of 
the near and flow field flapping on the far field pres- 
sure. Note in particular the alternating intervals in 
which the angular location of the larger values of 
p switches from xp = 0 to xp = n. The figure illus- 
trates that the timescale of the switching is generally 
greater than the timescales associated with the os- 
cillations in p as t varies. Furthermore, the figures 
illustrate the larger radiated pressure at mid angles 
as opposed to upstream or vertical angles. Finally, 
the figures illustrate the predominant low frequency 
beaming at mid angles (alternatively the preferred 
beaming of high frequencies upstream and at vertical 
angles). 

In Figure 6 we exhibit the power spectral den- 
sity (PSD) for p plotted against Strouhal number, 
St = /D/D e , where / is the frequency for a fixed 
location in the flow field. Specifically, we consider 
p at the point z/D = 1 . 1 , r/D = 0.5 and xp - 0 
which is in the shear layer of the jet. The PSD is 
taken over a window which excludes the initial pulse. 
The figure shows that most of the spectral energy is 
concentrated in a Strouhal number band between 
0.3 and 0.6. This is a higher range than the peak 
radiation observed in the far field, indicating that 
high frequency disturbances in the flow field are not 
necessarily efficient radiators of sound. Note that 
the PSD exhibits 2 distinct peaks at apparently in- 
commensurate frequencies, a result contrary to that 
found for 2 dimensional Cartesian jets with a sim- 
ilar Mach number where the flow field pressure ex- 
hibited narrow peaks at one particular frequency to- 
gether with harmonics. 6 We also note the similarity 
between Figure 6 and measured spectra in nonax- 
isymmetrically excited jets. 28 At this point in time 
we have not as yet conducted a systematic investiga- 
tion of the spectra at other points in the flow field. 
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4 Conclusion 

We have developed a code to solve the nonlinear 
Euler equations for nonaxisymmetric disturbances 
in an axisymmetric jet. The code is second order 
in time and fourth order in space and is written 
in cylindrical coordinates with the cylindrical axis 
the axis of the jet. The simulation includes both 
the near and far field of the jet as well as the up- 
stream region including a straight pipe from which 
the jet exits. The jet is excited by a transient, non- 
axsymmetric ring source located near the nozzle exit. 
The long time behavior of the jet is dominated by a 
train of disturbances generated at the nozzle lip and 
propagating downstream along the jet. These dis- 
turbances in turn generate sound as they grow and 
decay in the jet shear layer. 

The primary nonaxisymmetric feature that we 
have studied is a flapping mode of propagation 
within the jet. In these modes regions of high pres- 
sure switch 180° at certain time intervals. These in- 
tervals are generally long compared to the timescales 
of oscillation of the pressure within and outside of 
the jet. The flapping is also apparent in the pressure 
outside of the jet. 

There is a significant beaming of low frequency 
sound at mid angles to the jet axis, leading to an 
overall peak in the sound pressure level. In con- 
trast, high frequencies are beamed predominantly at 
angles of 90° and greater to the jet axis. 

The unsteady pressure in the jet flow field is 
primarily concentrated within the Strouhal number 
band between 0.3 and 0.6. Within this band there 
are 2 incommensurate spectral peaks. This is in 
contrast to results for 2 dimensional Cartesian jets 
where narrow peaks at a particular frequency and 
its harmonics are observed. 
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Pressui e contours bi sue 








Figure 3. Visualization of pressure in jet cross section for times t a - tf. Axial location 
z/D = 4.89. Radius of each cross section is approximately 2D. 








